Scaling in the one-dimensional Anderson localization problem in the region of 

fluctuation states. 
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We numerically study the distribution function of the conductivity (transmission) in the one- 
dimensional tight-binding Anderson model in the region of fluctuation states. We show that while 
single parameter scaling in this region is not valid, the distribution can still be described within 
a scaling approach based upon the ratio of two fundamental quantities, the localization length, 
hoc, and a new length. Is, related to the integral density of states. In an intermediate interval of 
the system's length L, lioc <^ L <^ la, the variance of the Lyapunov exponent does not follow the 
predictions of the central limit theorem, and may even grow with L. 



Introduction. A coherent transport in mesoscopic dis- 
ordered systems is characterized by strong fluctuations 
and non-self-averaging nature of the transport coeffi- 
cients such as conductance, g, or transmittance, T, [llQ- 
Therefore, a description of the transport in such systems 
requires dealing with entire distribution functions of the 
respective quantities. The scaling approach to the trans- 
port allows one to introduce a reduced "macroscopic" 
description of such distributions independent of micro- 
scopic details of the underlying Hamiltonians 0, 0] with 
parameters of the distributions playing the role of the 
scaling variables 0] . If the entire distribution can be pa- 
rameterized by a single parameter, the respective system 
is said to obey single parameter scaling (SPS). A quan- 
tity, which is most convenient to work with when describ- 
ing the statistics of transport is the Lyapunov exponent 
(LE), 7(L) = (l/2L)ln(l + 1/g), where L is the length 
of the system Finite size LE, 7(i), is self-averaging 

(approaches a non-random limit 7, when L ^ 00 0]), and 
its distribution approaches a Gaussian form for asymp- 
totically long systems. The Gaussian distributions are 
characterized by two parameters: the mean value, 7, and 
the variance, tr^, and the SPS hypothesis suggests that 
they are related to each other in a universal way. Such a 
relation, which can be expressed in the form 



1, 



(1) 



where hoc — I/7 is identified with the localization length, 
was first conjectured by Anderson et al. and repro- 
duced later by many authors within the framework of the 
phase randomization hypothesis 2] . The phase random- 
ization was proven rigorously for in-band states (those 
belonging to the spectrum of underlying ordered sys- 
tems) for some one-dimensional models (Anderson model 
and a continuous model with a white-noise random 
potential 0) as well as for some quasi-one-dimensional 
models 0. For the Lloyd model, the authors of Ref. 9 
showed that Eq. ^ (corrected by the factor of two) holds 
for in-band states even though the distribution of phases 
is not uniform. 

At the same time, numerical results presented in 
Ref. showed that Eq. is not valid for fluctuation 



states arising due to disorder outside of the original spec- 
trum. A boundary between SPS and non-SPS spectral 
regions in the exactly solvable Lloyd model was shown 
to be determined by a relation kodE) ~ h(,E) |^, where 
Is is a new length, defined through the number of states, 
N{E), per unit length, between E and the closest genuine 
spectral boundary: 



i: 



sin [7tN{E)] 



(2) 



In the region of fluctuation states, when N{E) <C 1 or 
1 - N{E) < 1, and Is > hoc, SPS is not vahd. Compli- 
menting analytical calculations by numerical simulations, 
the authors of Ref. 9 showed that the criterion for SPS 
found for the Lloyd model is valid for other models as 
well. 

Thus, it is clear that the problem of the statistics of 
transport in the region of fluctuation states requires a 
separate consideration. The distinction between this sit- 
uation and the case of in-band states can be qualified as a 
difference between under barrier tunnelling and over bar- 
rier scattering. In 3D this difference is clear: the latter 
case corresponds to the spectral region of extended states 
with the diffusive transport, while the former takes place 
in the region of localized states. The problem of the 
under barrier tunnelling in disordered systems was first 
considered in (for reviews of subsequent papers see 
Refs. In ID situation all states are localized, 

and the transmission for all energies can be described 
as a resonant tunnelling via rare transparent configu- 
rations [l3|- Therefore, the difference between the two 
transport regimes is more subtle and was noticed only 
recently 0, 0| . Correspondingly, while the case of the 
in-band states can be considered settled by the SPS the- 
ory, the properties of the distribution function of con- 
ductance/transmittance for pure one-dimensional case of 
under barrier tunnelling are studied very little. Besides 
obvious fundamental importance, an additional motiva- 
tion to deal with this problem comes from the devel- 
opment of photonic band gap materials, in which new 
type of fluctuation photonic states is possible ■ These 
states form "Lifshits tails" in the band-gaps of disordered 
photonic structures, and provide a unique opportunity to 
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study resonant under barrier tunnelling with scattering 
of light. 

The main objective of the present paper is to study 
numerically the distribution function of LE in the region 
of fluctuation states using ideas of the scaling approach. 
The main question, which we seek to answer is the follow- 
ing: "Is the distribution function of the LE in the non- 
SPS region determined completely by microscopic details 
of the respective Hamiltonian, or can it still be described 
macroscopically in an universal manner?" We show that 
the distribution of conductance in this region, while not 
completely universal, still demonstrates surprising scal- 
ing properties. In particular, using Monte Carlo sim- 
ulations for the one-dimensional tight-binding Anderson 
model with diagonal disorder, we find that for sufficiently 
long systems the function r, introduced in Eq. J^l, de- 
pends upon a single parameter, k = hoc/h- We also 
find strong deviations of the distribution function from 
the Gaussian form. However, the third moment, turns 
out to have the same scaling behavior as the variance, 
indicating that despite the deviation of the distribution 
function from the Gaussian, it still can be parameterized 
by only two parameters, koc, and Ig- 

Model and technical details. We consider the tight- 
binding model with a diagonal disorder, which is de- 
scribed by the following equations of motion 



(3) 



where random on-site energies C/„ are described by a 
uniform probability distribution: P(C/„) = 1/{2U) if 
\Un\ < U, and P{Un) = otherwise. LE is defined as 



-l{E) = lim - log llTjv • • • Till = hm ^{L), (4) 

N^oo 1\ L^oo 

where Th are transfer matrices 



Tk = 



E-Uk -1 



1 







(5) 



LE is calculated iteratively using Eq. in a standard 
way 0. To investigate the statistics of 7(L) in sys- 
tems with the finite length L, we keep the length finite 
and fixed while collecting statistics from about 120,000 
realizations. The integral density of states for each real- 
ization was calculated using the phase formalism Q and 
was averaged over all realizations. The resulting value 
was used to calculate the length Ig according to Eq. 
Studying the dependence of the distribution of 7 on L 
we take care to have L ^ koc for all strengths of the dis- 
ordered potential, and values of energy, E. However, 
in the region of fiuctuation states, where hoc < Is, it is 
possible to have Is > L ^ lioc- In this regime, which does 
not exist in the SPS region, the L-dependence of the vari- 
ance may be different from standard behavior given by 
the central limit theorem. In order to verify this assump- 
tion, we considered systems with lengths satisfying both 
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FIG. 1: Dependence of the scaling parameter t on k — hoc/h 
for a set of potentials {U — 0.08 - 0.155). On the insert 
dependence t(k) in the non-SPS region (k <C 1) is shown in 
log-log scale for 1 < L/ls < 5. 

L < Is, and L > Ig. When collecting statistics, we dis- 
carded all data corresponding to koc < 5, and Is > lOOOL. 
This way we ensured that our results are not influenced 
by non-representative fluctuations, and states localized 
over microscopical regions of the sample. 

Results. Fig. ^ shows the dependence of the scaling 
fimction t defined in Eq. on k in the asymptotic limit 
of very long systems, L ^ Is, hoc- Data used to gener- 
ate this figure were obtained for different values of E 
and U, and one can see that they all fall nicely on the 
scaling curve t(k). For k > 1, t approaches its uni- 
versal SPS value of unity while for smaller k it steeply 
decreases. A similar result was obtained for a periodic- 
on-average model in Ref. 9] , where the scaling function 
t(k) was originally proposed. Our results convincingly 
show that (7^ can indeed be expressed in terms of the 
scaling function t(k) regardless of the microscopic na- 
ture of the model under consideration. 

While the exact shape of the function t{k) varies 
slightly from model to model (compare to results of 
Ref. the essential qualitative properties of r seem to 
be quite universal. We are most interested in the region 
K < 1, where r demonstrates a sharp decrease. Ana- 
lytical calculations carried out in Ref. for the Lloyd 
model produced t(k) w k. Our results show, however, a 
much steeper decrease of r. Also, in the model consid- 
ered here, t{k) must remain non-zero for k = 0. Indeed, 
K = corresponds to the exact genuine boundary of the 
spectrum of our system. Unlike the Lloyd model, where 
the spectrum boundary is at infinity, in our model the 
boundaries of the spectrum are at Ei, = ±{2 + U). The 
variance of LE does not vanish at finite energies, and 
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FIG. 2: Index of the scaling parameter a (filled squares, left 
axis) and factor C (circles, right axis) as functions of L/h. 

therefore = t(0) is not equal to zero. 

In order to understand the behavior of t at k ^ 1 
we conduct a detailed study of this region for systems 
with different L. Our results can be summarized in the 
following form 



= Ck" 



(6) 



Replotting t{k) in the log — log coordinates for k < 1 (see 
insert in Fig.^l we see that while k changes by more than 
two orders of magnitude, the data form a good straight 
line with the exception of points corresponding to ex- 
tremely small values of k. This means that can be 
neglected for the most of the non-SPS region, and be- 
comes significant only in the immediate vicinity of Eh- 
Using linear regression we can estimate parameters C 
and a for systems with different lengths. The results of 
the fit reveal that C and a are constants independent 
of any parameters of the system under consideration for 
L > Is- This result confirms the one-parameter form of 
t(k) given by the first term of Eq. (j^ for sufficiently long 
systems. The degree of universality of these coefficients 
still remains an open question requiring similar studies 
of other models. We can speculate, however, that it is 
likely that systems can be divided into several universal- 
ity classes on the basis of the values of C and a. 

For shorter systems with L < Is both C and a show ad- 
ditional dependence upon the length L, see Fig. |21 where 
a and C are plotted versus L/lg. 

The results of this analysis lead to two important con- 
clusions. First, the length Is not only establishes the 
boundary between SPS and non-SPS regions of the spec- 
trum, but also determines a crossover system length 
marking the transition to systems with a universal single 
parameter form for t. Second, as it was anticipated, in 



the regime koc L < Is the scaling of the variance of LE 
changes from the simple 1/L dependence to more a com- 
plicated form due to the dependence of a on L/ls- We 
attempted to fit this dependence in the region of small 
L/ls by several types of trial functions; the best fit was 
obtained with a{L/ls) ^ In {Is/L). With this assumption 
a new scaling for the variance becomes 



cx 



1 



■exp[a(L//,)lnK] cx L-(i+'"'"'. 



(7) 



It is interesting to note that when k decreases, 1 + In k 
becomes negative and cr^ starts growing with L in this 
interval of lengths. This behavior can be qualitatively 
understood from the following arguments: The condi- 
tion L <ti Is means that for the most of the realizations 
of the random potential no states exist in the energy in- 
terval under discussion. The transmission through such 
realizations fluctuates rather weakly. The greatest con- 
tribution to the transmission fluctuations give those few 
realizations that can support at least a single state. The 
probability for such realizations to arise, grows when the 
length of the system increases, resulting in the respective 
increase of a^. 

This behavior, of course, breaks down for very large 
values of Is , which correspond to states close to the gen- 
uine spectral boundary, because for these states is 
determined by a non-universal correction to r given by 
Tiim- This limiting value can be found using the weak 
disorder expansion of Derrida et al. [l^ , which gives 



Tli' 



47osinh 70 



cx VU, 



(8) 



where afj — is the variance of the potential, 70 

is LE in the gap region of the original ordered system, 
and we use the fact that Ei, = 2 + U for the system 
under consideration. One can see from this expression 
that Tiim depends on microscopic characteristics of the 
original Hamiltonian. However, for a model with the 
Gaussian distribution of site energies, the genuine spec- 
tral boundary lies at infinity, where 70 — > co. The first of 
the equalities in Eq. (jS)) , which can be applied to various 
distributions, in this case gives Tum — 0. We can expect 
this to be true for all models with spectral boundaries at 
infinity. For this class of models, t{k) gives a completely 
universal, at least within a given model, description of 
the variance of LE. 

In order to characterize deviations of the distribution 
function of LE from the Gaussian form, we also studied 
scaling properties of the third cumulant g = ((7— (7))'^), 
which describes the skewness of the distribution func- 
tion. To analyze scaling properties of g wc introduced a 
function analogous to r 

T3 = gL^koc- (9) 

One can see from Fig.|3|that while data for the param- 
eter Ta are rather noisy, it shows a relatively good scaling 
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FIG. 3: Dependence of ra on kappa deeply inside of the non- 
SPS region (k <C 1). On the insert scaled third cumulant gL^ 
is depicted as a function of the energy near the band edge for 
the potential U = 0.05 and sample lengths from 1040 to 2160. 



behavior as a function of the single parameter k, in the 
non-SPS region. This fact itself is quite remarkable since 
it demonstrates that even deviations from the Gaussian 
in the region of fluctuation states can be described within 
the scaling procedure suggested here. It is interesting to 
note that the sign of the skewness changes not very far 
from the boundary between SPS and non-SPS regions 
of the spectrum. In the SPS region the absolute value 
of skewness decreases dramatically becoming essentially 
zero within the accuracy of our calculations. This is illus- 
trated in the insert in Fig. 13 which represents the energy 
dependence of the third moment of the distribution. The 
skewness also decreases deeper in non-SPS regions where 
it becomes extremely small beyond the genuine spectral 
boundary. The quality of our raw data did not allow us 
to determine whether also depends on L/ls for shorter 
systems, but we expect that T3 behaves similar to r in 
this regard. 

Conclusion. In this paper we carried out a detail 
study of the distribution function of conductance in the 
spectral region of fluctuation states. We showed that 
apart from a small non-universal contribution, which is 
only important in the immediate vicinity of the genuine 
spectral boundary, the conductance distribution in this 
region can be described using a simple scaling approach. 
Within this approach the variance is described by the 
scaling function r, Eq. (Q, and the third moment of the 
distribution is characterized by the function T3, Eq. I^. 
For long enough systems both scaling functions depend 
on the single variable k — hoc/h- The presence of such a 
scaling behavior would be natural for the model with the 



white noise potential, because such a model has a natural 
scaling variable E/u^'^ 0, and our scaling parameter k 
depends upon this only variable Our numerical re- 
sults showed, however, that the parameter k provides a 
more universal description of the distribution function 
valid also outside of the white-noise model. While we 
only considered here the tight-binding model, we believe 
that our results qualitatively describe statistics of light 
transmission through band-gaps of disordered photonic 
crystals. Experimental measurements of the transmission 
distribution in such systems can be used for verification 
of our results and as a method of measuring parameters 
hoc and Is- 
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